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The general equations of motion for two dimensional singular Laplacian growth are derived using 
the conformal mapping method. In the singular case, all singularities of the conformal map are on the 
unit circle, and the map is a degenerate Schwarz-Christoffel map. The equations of motion describe 
the motions of these singularities. Despite the typical fractal-like outcomes of Laplacian growth 
processes, the equations of motion are shown to be not particularly sensitive to initial conditions. 
It is argued that the sensitivity of this system derives from a novel cause, the non-uniqueness of 
solutions to the differential system. By a mechanism of singularity creation, every solution can 
become more complex, even in the absence of noise, without violating the growth law. These 
processes are permitted, but are not required, meaning the equation of motion does not determine 
the motion, even in the small. 



I. INTRODUCTION 



Laplacian growth - growth of a region D along the 
gradient of its external Green's function - is a model 
for a number of growth processes which occur in nature, 
among them growth by electrodeposition diffusion- 
limited aggregation j2|, and viscous fingering at fluid- 
fluid interfaces [||. These natural processes exhibit 
very complicated morphologies, as does the mathemati- 
cal model. In spite of the large amount of work which 
has been done, one still has the feeling that there is some- 
thing mysterious about Laplacian growth. In particular, 
its extreme sensitivity to perturbations makes it difficult 
to interpret experiments, real or numerical. 

A special role is played in the two-dimensional prob- 
lem by the conformal mapping method, as developed in 
0, 0, 1], @, @, for example. This method gives in- 
sight difficult to attain otherwise. It is especially sim- 
ple in radial geometry In outline, one parametrizes the 
growing 2-dimcnsional region D, thought of as occupy- 
ing a bounded, simply connected region in the complex 
w-plane, by the conformal map 

w = G{z) (1) 

which takes the exterior of the unit disk, \z\ > 1, onto the 
exterior of D. The growth of D is then represented by 
the time dependence of the conformal map G. Since G 
is holomorphic in \z\ > 1, and its dependence at infinity 
is also prescribed, G may, in turn, be parametrized by 
its singularities in the unit disk \z\ < 1. The growth be- 
comes the dynamics of those singularities. This method 
eliminates sources of inaccuracy which are unavoidable 
in other methods, for example the statistical noise which 
accompanies DLA simulations, or some of the roundoff 
and truncation errors of more straightforward integration 
methods. That does not make the conformal mapping 
method necessarily more realistic, of course. Indeed, the 



noise in other numerical methods may model actual phys- 
ical noise, and hence be desirable, if one's aim is to model 
particular examples of Laplacian-like growth. With the 
conformal mapping method we aim rather to strip the 
problem down to its simplest form, to see what remains 
and is common to all such processes. 

In reference this process of stripping down was taken 
one level further, and an unexpected phenomenon came 
to light. This case might be called singular Laplacian 
growth, because it is the limiting case in which all the 
singularities of G are on the unit circle \z\ — 1. In this 
case G degenerates to a Schwarz-Cristoffel map onto the 
exterior of a degenerate polygon D of zero area (i.e., D 
looks like a tree graph). In this limit the dynamics be- 
comes 1-dimensional, and can be understood completely 
The surprise is that the dynamics allows the singularities 
of G to split and proliferate, but it does not require this. 
That is, while the dynamics is formally given by differ- 
ential equations, the solutions, for given initial data, are 
not unique. The comment was made in Ref. that this 
appears "unlike any other physical model." In particular, 
it is not the same as being very sensitive to initial con- 
ditions, as one might have assumed. In fact, as we show 
below, singular Laplacian growth is not at all sensitive in 
this way. Its peculiarities have a different origin. 

Reference [|7| gave only the simplest example (in which 
all maps could be written down explicitly), and not a 
general computable theory. The present paper gives the 
general theory. 

II. DYNAMICS OF SINGULARITIES 

Let w = G(t, z) be a time-dependent conformal map 
of the exterior of the unit disk in the z plane onto the 
exterior of the domain D in the w plane. G has the form 
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G{t,z) = zJ2c k (t); 
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Suppose for the moment that OD, the boundary of D, is 
an analytic Jordan curve so that there is no difficulty in 
defining 



g(t,6)=G(t,e l9 ). 



(3) 



As shown by Shraiman and Bensimon j5|, Laplacian 
growth implies that the conformal map has time depen- 
dence given by 
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As a conformal map, H has singularities at points on 
the unit circle which we have called f3j and 7^. The 
image of an arc of the unit circle under H, so long as 
it does not contain a singularity, is a straight line seg- 
ment. At the singularity (3j, however, the image line 
turns through the angle ajw, which may be either posi- 
tive (counterclockwise) or negative (clockwise). It is un- 
derstood that \ctj\ < 1. At the singularity jk the image 
turns through the angle n, i.e., the line retraces itself: see 
Fig. 1. The 7 singularity may seem to be only a special 
case of the (3 singularity, corresponding to a = 1, but we 
have distinguished it because growth takes place entirely 
at the 7 singularities: the 7's are different. (Notation: 
a= "angle"; (3= "branch point" ; 7= "growth tip." The 
ctj of this paper was called a 3 - — 1 in 0.) 

The dj, /3j, and jk are by no means arbitrary. First, 
because the image polygon turns through a total angle 
2ir, by conformality, one must have, since there are M 
branchpoints and N growth tips, 
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Define new scaled variables 

afc = Cfc/co, b k = d k /d 
and a new independent variable s(t) by 

ds/dt = do. 
In terms of these variables, Eq. (^) becomes 



^ = (fc-2K + 2^(l-j)a,& fc _,, 
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Because of the rescaling in Eqs.(^) and @, Eq. (||) con- 
tinues to make sense even in the limit as singularities of 
the conformal map move onto the unit circle. For ex- 
ample, \bk\ < 1 for all k, even though dk blows up. We 
define the bkS, in case there are singularities on the unit 
circle, to have their limiting values as the singularities 
move onto the unit circle from the inside. With this un- 
derstanding, Eq. (||) describes Laplacian growth, both 
singular and nonsingular, in terms of the scaled mapping 
function 



Second, because the image polygon looks like a tree 
graph, each edge of the image is traversed twice, once in 
each direction. This means that the integrals of OH/Oz 
along arcs of the unit circle from one singularity to the 
next, which are singular integrals, must cancel in pairs, 
an intricate condition on the positions of the /3's and 7's. 

Suppose at time s = we have such a conformal map 
H. The equation of motion for H, Eq. (^), should be 
recast as equations of motion for the singularities of H. 
The cifc's of Eq. @ are just the coefficients in the power 
series for H, according to Eq. ([To|), but we still need 
the 6fc's. To compute the 6^'s, using Eqs. (JsJ) and (0), 
we must Fourier transform \0H/0z\~ 2 , restricted to the 
unit circle, with the singularities displaced slightly inside. 
The Fourier transform integrals are dominated by the 7 
singularities, and as the singularities move onto the unit 
circle, the entire contribution comes from them. Thus 
there is a simple formula for bk, 
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In the singular case, in which all singularities of H are 
on the unit circle, H is a Schwarz-Christoffel map onto a 
degenerate polygon, and therefore its derivative has the 
form 



Here the "weights" Vj are determined by 
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with the constant of proportionality determined by the 
normalization 
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This is gradient flow in the space of parameters 

(15) 

(/3i, /3m, 7i, 7tv) endowed with the metric tensor 



Define the function 
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Then multiplying each side of Eq. @ by 2 k+1 and sum- 
ming over k gives 
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It is dH/dz rather than H that we know explicitly, so 
take the derivative of Eq. (|l7|) with respect to z. The 
left side becomes 
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and the right side becomes an explicitly known expres- 
sion. One can now cancel the factor —idH/dz in all 
terms. Multiplying by (1 — e l/3j / z) and taking the limit 
as 2 — > e ll3j isolates df3j/ds, and similarly multiplying 
by (1 — e m I z) and taking the limit as 2 — > e llk isolates 
drfk/ds. The result, after algebraic simplification, using 
Eqs. @ and ©, is 
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We can also note 
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It is understood in Eq. ( |22| ) that is real and posi- 
tive. Eqs. (|l9|)-(p3|) represent the dynamics of singular 
Laplacian growth as an autonomous system of ODE's. 
Remarkably, this system is a kind of gradient system: 
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g = diag(ai, ...,a M , 1, 1) 



(26) 



Since, according to Eq. (|12|), the a's are negative, on 
average, if N > 2, this metric is indefinite. The gradient 
flow is toward certain critical points of in W which are 
not minima. These critical points are the equilibria of 
singular Laplacian growth. They can be found by inte- 
grating the system of Eqs. (|l9[)-(|23l). Even if the starting 
state does not satisfy all the conditions described after 
Eq. (|l2|), it will still approach a state which does sat- 
isfy them. We describe these equilibrium states more 
precisely below. This stability of the flow, which is a 
familiar property of gradient flows, is an indication that 
singular Laplacian growth is not sensitive to initial con- 
ditions, contrary to what one might have expected, and 
hence that the peculiar sensitivity of Laplacian growth 
in general has either been lost in the passage to the sin- 
gular case, or that it arises from some other cause. We 
suggest below that that other cause is the non-uniqueness 
property of the system, still to be described. 

The dynamics of singularities described by Eqs. (|l^)- 
( p3| ) can be pictured very simply. The /3's are repelled by 
the 7's on the unit circle, and the 7's repel each other. 
The "strength" with which each 7^. repels other singular- 
ities is given by the corresponding v k (always > 0). The 
/3's, on the other hand, do not interact directly with each 
other. No singularity can pass through another one - 
they always keep the same order around the unit circle. 
Tljipj^ /3's between any two adjacent 7's are, however, 
driven by them toward some intermediate point where, 
in effect, they coalesce into a single effective (3, charac- 
terized by a single effective a which is the sum of all the 
contributing a's. The way a (3 singularity approaches its 
limit position is the way x(s) approaches in 



namely 
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x(s) = 2:0/(1 + x s), 
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that is, it takes infinite time. By the second derivative 
test, there is only one equilibrium position for the /3's be- 
tween each pair of 7's. Here all the /3's will collect. Thus 
in the limit as s — > 00 the equilibrium states of singular 
Laplacian growth have 7's and f3's alternating, and look 
like N needles radiating from a single central point. One 
can even write a formula in closed form for H in this case, 
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where the /3's and a's are the effective ones. We can 
also understand this outcome, in a more physical way, by 
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realizing that the continual rescaling means all internal 
structure shrinks to a point, leaving only the growth tips 
as visible features. What is not obvious from this descrip- 
tion, but is observed, is that typically some of the 7's are 
entrained with the /3's, and coalesce with them (where 
they contribute +1 to the effective a). This amounts to 
the scaling away of needles. It turns out that the generic 
stable equilibria have N < 3. If the initial configuration 
has N > 3, and is the least bit asymmetrical, some of 
the growing tips lose out in the competition to grow, and 
disappear as s — > 00, leaving only three (or fewer) needles 
in the limit. 

This result might appear puzzling, since it seems to 
imply that Laplacian growth should be a process of sim- 
plification, contrary to the increasing complexity which is 
observed, and which is the whole motivation for studying 
it. That puzzle is resolved in the next section. 

III. NON-UNIQUENESS OF THE DYNAMICS 

If one reverses the sense of time and integrates the 
system backwards, the repelling character of the 7's be- 
comes an attraction. In particular, two /3's adjacent on 
either side of a 7 may be attracted to it, move toward 
it, and coalesce with it, essentially annihilating, leaving 
just the 7. Unlike the coalescence described at the end of 
the previous section, which takes infinite time, this coa- 
lescence occurs in finite time. Roughly, one can estimate 
from Eq. (|l9|) that a (3 approaches a nearby 7 the way 
x(s) approaches in 

dx/ds = 1/x, (30) 

(integrating backward from s = 0), namely 

x(s) = yja* + 2a. (31) 

At times earlier than the coalescence time — Xg/2, the /3 
singularities were simply not present. If one now exam- 
ines this solution to the system with the usual forward 
sense of time, one sees, at some arbitrary time -Xq/2, two 
j3 singularities suddenly produced on either side of a 7, 
which hadn't been there before. To satisfy Eq. (fl2|), the 
ct's which characterize these /3's must add to zero. The 
geometrical effect of this process is that a kink of de- 
viation angle a suddenly appears in the growing needle 
represented by 7, like the kink shown in Fig. 1, which 
might have formed from a single straight needle. This 
kinking may happen at any arbitrary time. A more care- 
ful argument (in the Appendix) says that if a kink forms 
at 7 at s = 0, the leading behavior in the motion of 
singularities is 

V 1 + OL2 

(7-/3 2 )cc-,/i±!V/ 2 (33) 
V 1 + CX\ 



with a± — —ci2- 

In addition, a second kind of coalescence is seen in 
backward integration, in which two /3's, with angles ct\ 
and CH2 on either side of a 7, coalesce to leave a single /3 
with angle l+a^+o^. This happens only if l+ai+a 2 > 
and ai + C12 < 0. Geometrically it corresponds to the 
shrinking away of a needle in finite time (the growth tip 
7 is lost), on the outside of a corner of angle 1 + a\ + 0.2- 
What it means in forward integration is that at any time 
a needle may begin growing on the outside of a corner, as 
in the process which takes Fig. 1 to Fig. 2. The motion 
of singularities in this case, in leading order, is 



7 _^ oc ./i±^l s V2(l+« 1 +a 2 ) (34) 

V 1 + Oil 

j-fooc- /l±^I s i/2(i+« 1 +« 2 ) (35) 

V 1 + Oil 

with 1 + oci + CY2 > 0, cxi + 0L2 < 0. (In Ref. [Q] the factor 
1 + ai + 012 in the exponent's denominator mistakenly 
appeared in the numerator.) 



IV. DISCUSSION 

The observations of Section III mean that the system of 
ODE's Eqs. although it looks unremarkable , 

has the peculiar property that its solutions are highly 
non-unique. New singularities can appear by the above 
two elementary processes at any time. In combination 
one has more complicated processes: a kink followed by 
a new needle at the outside of the new corner amounts 
to tip splitting, for example, and this can happen at any 
time. The equilibria described in Section II are never 
attained if such processes, which are allowed by the dif- 
ferential equation, continually intervene. Thus singular 
Laplacian growth supports complex non-equilibrium be- 
havior after all. 

It is interesting to see what the model looks like if one 
integrates it forward, introduces new singularities, inte- 
grates again, adds more singularities, etc. Examples are 
shown in Figs. (3) and (4), where symmetrical tip split- 
ting was introduced at intervals of 0.1 time unit. To inter- 
pret the evolving positions of singularities in terms of the 
corresponding image region D, which is what is shown, 
it was necessary to integrate Eq. ( pd| ) numerically. Each 
edge is represented by a singular integral. These integrals 
were done by Gauss-Jacobi integration, as described by 
L.N. Trefethen in Ref. The accumulating error in 
these numerical integrals, as one steps along each edge 
to the next, especially in light of the usual sensitivity of 
numerical conformal maps, might have been expected to 
produce nonsensical pictures, but in fact the numerical 
error (failure to retrace edges accurately) is just barely 
visible in these examples. (Eventually, of course, the ac- 
cumulating error does become large, but the good numer- 
ical behavior of the system again makes the point that 
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singular Laplacian growth is not particularly sensitive to 
error or noise. Its sensitivity to perturbations comes en- 
tirely through the non-uniqueness property.) 



V. RELATION TO OTHER WORK 

Most of those who have used the conformal mapping 
method have followed Shraiman and Bensimon || in re- 
stricting the derivative of the conformal map H to be 
a rational function. From some points of view this is a 
rather drastic restriction on the analytic structure of H. 
Whether it is a good enough representation of H to learn 
the full implications of the conformal mapping method is 
not clear. Arguments that the boundary value of H' can 
always be approximated by the boundary value of a ra- 
tional function are not very convincing in a context where 
it is precisely the nature of the singularities which is the 
basis of the theory. It had already been noticed in Ref. |t]] 
that branch points play an essential role in the singular 
theory. Nonetheless, an interesting comparison between 
the singular case and the rational case is possible. 

An example is Ref. in which R. Blumenfeld and 
R.C. Ball invent a mechanism of "particle creation" (i.e 
singularity creation) to model tip splitting. In their 
model, since H' is rational, the only singularities are the 
zeros and poles of H' . The mechanism they propose is 
that a zero creates a second zero and a pole. The two 
zeros represent the two growing tips after the split, and 
the pole represents the division between them. 

Tip splitting in the singular theory, as described in Ref. 
and in this paper, does not have to be invented - it 
is naturally and unavoidably part of the theory - and it 
looks slightly more complicated: a 7 gives rise to three 
/3's and another 7 (as in Figs. 1 and 2). But in fact this 
amounts to the same thing. The resulting two 7's are ze- 
ros of H', and, by the geometry of the situation, the three 
new a's add to —1. This means that the three /?'s, from a 
distance much greater than their mutual separation, look 
like a pole of H' (see Eq. (11)). The mechanism proposed 
by Blumenfeld and Ball is thus a kind of smeared version 
of the singular mechanism, already described in Ref. 

It is especially remarkable that Blumenfeld and Ball 
invented their mechanism entirely on the basis of physi- 
cal phenomenology, and were unaware of Ref. . Their 
mechanism of particle creation, although it is ad hoc, cor- 
responds as precisely as it could have to the only mech- 
anism in the singular theory for non-trivial dynamics. 
This suggests that the singular theory is close enough to 
real phenomenology to be useful, and that it does retain 
the essential features of Laplacian growth. 



VI. GENERALITIES 

To focus on the details of singular Laplacian growth is, 
to some extent, to sidestep a much bigger question: what 



is going on here with non-uniqueness? Aren't differential 
equation supposed to have unique solutions? We all know 
textbook examples where uniqueness fails, but the fail- 
ure occurs on some small set, and for equations which 
wouldn't arise in physics. Here are equations which arise 
in a system which has been much discussed in physics, 
and uniqueness fails for every solution at every time. The 
least one can say is that the equations of motion do not 
determine the motion, even in the small. 

I believe this is actually mathematical terra incognita. 
Such equations do not even have a name. How would 
one characterize them generally? Are they in some sense 
common, or are they rare? I think of calling them "frag- 
ile differential equations," because, at least in this ex- 
ample, the non-uniqueness arises by the tendency of sin- 
gularities to "break apart," but perhaps a more general 
understanding would reveal that this name is somehow 
misleading. "Fragile" sounds a little bit like "fractal," 
but is not the same, another reason I like the name. 

On a more physical level, what does it mean for a phys- 
ical system if it is described by equations which, in some 
limit, become "fragile"? A fragile system does not fully 
determine the evolution, but it does restrict it. What 
is the nature of the restriction? These seem like good 
questions for the future. 



APPENDIX 



We derive Eqs. (|32j)-(|35|), the leading behavior of sin- 
gularities [3\, 7, /?2 when they are very close to each other 
(in that order), and not close to other singularities. Let 
«i and 012 be the corresponding angle parameters, and v 
the "strength" of 7. According to Eqs. (p^|)-(pO|), keep- 
ing only the most singular terms, in leading order they 
obey 
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P = 7 - 0i 
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Then, subtracting, 
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Dividing, we have the homogeneous equation 
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(37) 
(38) 



(39) 
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dP _ (l + a 1 )Q + a 2 P 
~d~Q ~ ai Q + (l + a 2 )P 



(43) 



which separates when written in terms of the variable 
P/Q. The complete solution, in implicit form, is 

(P - Qf 1 [(1 + a 2 )P + (1 + ai )Qp = const. (44) 

where 



v\ = (1 + ai + a 2 )/(2 + a.\ + a 2 ) 
v 2 = 1/(2 + «i + a- 2 ) 



(45) 
(46) 



Since this is supposed to hold as P, Q approach zero, the 
only relevant value of the constant is zero. The solution 
P = Q is not relevant to this situation, since P and Q 
must have opposite sign. Thus 



(l + a 2 )P+(l + ai )Q = 0. 



(47) 



Using Eqs. (J21|)- (p3|) together with the fact, found in Eq. 
([l?]), that 7 — f}\ and 7 — (3 2 are simply proportional, we 
see that v is non-singular if ct\ + a 2 > 0, and 



if ai + a 2 < 0. Thus from Eq. (|4 
dP/ds - P" 1 



(48) 



(49) 



if a>i + a 2 > 0, as in the rough argument of Section III, 
and 



^ p-2ai-2a 2 -l 

ds 



(50) 



if <x\ + a 2 < 0. These results are summarized in Eqs. 
©-©• 
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FIG. 1. The image under a degenerate Schwarz-Christoffel map of an arc containing singularities 7, /?2,--- The image 

comes in from the left, turns through an angle ain at H(/3i), stops and reverses direction at #(7), turns through an angle ain 
at H(P2), and exits on the left. In this example, showing a kink in a growing needle, ai = —02. 



H(y) 




H(y') 



FIG. 2. A new needle may grow on the outside of the corner in Fig. 1. Here the singularity f3\ has split into two branch 
points, /3[ and /V', and a new growth tip 7'. 
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FIG. 3. At each time interval As = 0.1 the growing tip with largest strength is split. 

random needles radiating from a point, but the growth law is completely deterministic. 
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The initial configuration was four 
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FIG. 4. At each time interval As = 0.1 a growing tip is randomly selected and split. The probability of a tip's being 
selected is proportional to its strength v. The initial configuration was three random needles radiating from a point. 
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